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Abstract 

We study the Z boson pair production mediated by the Kaluza-Klein (KK) graviton in large 
extra dimensions (LED) at the CERN Large Hadron Collider (LHC). We use the partial wave 
unitarity to discuss the constraints on the process energy scale in order to give a self-consistent 
calculation. We find that the LED contributions can enhance the Z boson pair production cross 
sections significantly when the fundamental scale Ms of the large extra dimensions is up to several 
TeV. We also show that the kinematic distributions of the LED signals are greatly different from 
the standard model ones and the LHC can probe the Mg values up to 4.3 ~ 5.6 TeV at 3a level 
depending on the number of the extra dimensions. 
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I. INTRODUCTION 



The idea that quantum gravity can appear at the TeV energy scale well below the Planck 



mass Mp ~ 1.2 x 10 19 GeV was proposed in the 1990's The large extra dimen- 

sions(LED) model [2j introduced by Arkani-Hamed, Dimopoulos and Davli has attracted 
much attention as the presence of large extra dimensions brings a new solution to the hier- 
archy problem. In the LED model, n extra spatial dimensions are compactified on a torus 
with common circumference R, and a 3-brane is introduced in which the standard model 
(SM) particles live. The SM gauge interactions are confined to this brane and only the 
gravity can propagate in the extra dimensions. Then the four-dimensional Planck scale Mp 
is no longer the relevant scale but is related to the fundamental scale Ms as follows J4]: 

Mp ~ M^ +2 R n , (1) 

where Ms ~ TeV. According to Eq. (CQ), deviations from the usual Newtonian gravitational 
force law can be expected at distances smaller than R ~ 2 x 10 -17+32 / n cm [4J. For n > 2, 
LED is consistent with the current experiments since gravitational forces currently are only 
well probed at distances about 40 /im [5] (However for n = 2, there are constraints arising 
from, e.g., supernova cooling {(|, which require Ms > 14TeV if n — 2). Recently, there are 
also some new constraints Ms > ITeV [7] from direct search at the Tevatron. 
The gravitational couplings of the SM fields are given by 



dxW-g£ SM (g,S,V,F), (2) 

where g is the induced metric in four dimensions. Expanding g with the gravitational fields, 
we obtain 

- 1 E / d ^Ch^%u + cu^n), (3) 

k 

where k = y/lQir/Mp, uo = a/2/ (3n + 6), T^ u is the energy-momentum tensor and h^' h , <fi k 
are the four dimensional Kaluza-Klein (KK) gravitons and KK scalars, respectively. From 
Eq. ([3]) it can be seen that the individual KK mode couples with a strength 1/Mp to the SM 
fields. However, since there are many KK modes, the total coupling strength is of the order 
I/M5 after summing up all of them. For the phenomenology studies, we use the convention 
introduced in Ref. 8|, where Ms is chosen to satisfy 

n 2 R n = 87r(47r) ri / 2 r(n/2)M- ( " +2) . (4) 



In this year, the LHC will start running and begin searching for the Higgs boson, further 
testing the SM, and also searching for new physics beyond the standard model. The Z boson 
pair production is an important process at the LHC and has been studied extensively. In 
the SM the Z boson pair can be produced directly through electroweak interactions with a 
total cross section of about 17 pb 9| at the LHC, corresponding to 1.7 x 10 6 events with an 
integrated luminosity of 100 fb _1 . With such a high production rate, it can serve as a probe 
for testing SM at TeV scale and also be an important background for new physics signals. 
Besides, the Z boson pair production can get additional contributions from new physics, 
for example, from the tree level processes in extra dimension models lOj, llll], little Higgs 



model 



12| and from loop level process in supersymmetry model 13j, thus it can be used as 



a probe of new physics beyond standard model. 

Usually, the most promising processes to test the LED model at the LHC are dilepton 
production, diphoton production, dijet production and single vector boson or jet production 



associated with KK graviton lQj, |14|]. However, the Z boson pair production also can be an 



important process especially for low Ms value, the total cross sections and invariant mass 
distributions of this process have been given in Ref. [lol |. In this paper we will further study 
the Z boson pair production through the s-channel KK graviton exchange at the LHC. We 
focus on the 4 charged lepton final states and perform a Monte Carlo simulation to see 
the new physics effects on several kinematic distributions. Finally, we also show the LHC 
discovery potential of M$ with different values of n. 

The arrangement of this paper is as follows. In Sec. II, we give the helicity amplitudes 
and discuss the unitarity constraints on the process energy scale. Sec. Ill shows the results 
of the total cross sections and Sec. IV contains the simulation results of several kinematic 
distributions and also the LHC reach of Ms- Sec. V is a brief conclusion. 



II. HELICITY AMPLITUDES AND UNITARITY CONSTRAINTS 

In the LED model the Z boson pair can be produced through KK graviton exchange from 
the gg fusion and qq annihilation, respectively, at the LHC. In order to calculate the relevant 
amplitudes we need to sum over all the KK graviton propagators, 

d( S ) = y — V— • ( 5 ) 

^ s - mt + ie 

k k 
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Since the KK graviton mass separation of 0(1/ R) is much smaller than all the other physical 
scales involved, we can obtain the above sum in the continuum limit Q], 

„n/2— 1 pn 

J(s) = r(,,./2)(4,)^ + aI < A/ ^ < 6 > 

where A ~ O(Ms) is the explicit ultraviolet cutoff scale of the effective field theory, and 

( n_1 )/ 2 1 / a \ 2fc-l , / . z - \ 

At the LHC the subprocess energy is comparable to M s , so we need to keep all the terms 
in Eq. ©. 

Since the scattering amplitudes grow quickly with y/s, we should cut off the subprocess 
energy at some scale below A in order not to violate the unitarity. And we use the partial 
wave unitarity to discuss the constraint on y/s. The J-partial wave amplitudes are given by 



647r 



J dcos9 d J ^,(cos9) [M XlX2X '" Xi ] , (8) 



where \i = Ai — A2, n' = A3 — A4, 9 is the scattering angle and d J t is the Wigner function [15 ]. 
For the subprocess gg — > ZZ, the individual non-vanishing helicity amplitudes are 



iM+-+ = 


-iM + - - 


= 2v / 2i5 crf M s " n " 2 m z 7rs :i * i (l + cos(0)) sm{6)A; 


iM+— = 


-iM+- Q+ 


= 2v / 2i5 crf M s " n ~ 2 m z 7rs^(l - cos(0)) sm(6)A; 


iM+-++ = 


iM+— = 


AiS cd Ms~ n ' 2 ra z 2 7is n/2 sin 2 (fl)„4; 


iM+-+- = 


At5 cd Ms- n - 


- 2 7r S t+i CO s 4 (9/2) A; 


iM+—+ = 




- 2 7rst +1 sin 4 (9/2) A; 


iM +-oo = 




- 2 7rs n/2 (4m z 2 + s) sin 2 (9)A, 



(9) 

where c, d are the gluon color indices and 

A = -2I(A/V«) + iTT. (10) 

Since the energy scale we considered here is much larger than mz, so we can neglect all 
the terms that are proportional to mz- Then the main contributions arise from A4 + h _, 



.M-i |- and .M+_ ocb an d the only non- vanishing partial wave amplitudes correspond to the 

J = 2 partial wave, which are given by 

i / r \ n+2 

a 2±2 = ^f^J A a 2 w = -a 2 22 /V§. (11) 

All the above amplitudes a| 2 , a|-2 an d a io contribute to the imaginary part of gg elastic 
scattering amplitudes according to the optical theorem, so the partial wave unitarity leads 
to 

kl 2 | 2 + |al- 2 | 2 + l«lol 2 < 1- (12) 

The partial wave amplitudes of qq — ► ZZ lead to weaker constraints compared with the 
ones of gg — > ZZ, and are thus not shown here. In Fig. [TJwe show the unitarity constraints 
on y^s with different A/M$ values according to Eq. ( fl2l) . We can see that the constraints 
on -y/i become strong for large A/ Ms values. In the following calculations we set A = Ms 
and choose the energy cut-off to be 0.9Ms for simplicity, which is allowed by the unitarity 
constraints. 



III. CROSS SECTION CALCULATIONS 



The additional contributions to the partonic differential cross sections of the Z boson pair 
production in the LED model are given by 

da qq - sy/1 - Ax a 2 (l + b 2 )T{n)T (n) 
dcosO ~ 24 

z 4 (1 - 4x) 2 + 10z 2 x (1 - Ax) - 2x (1 - Ax) - 1 
(1 - z 2 ) (1 -Ax) +Ax 2 

ttsVI -AxT 2 {n) (tt 2 + AT 2 (n)) 



and 



with 



384 

\2 i 2 



{ -3z 4 (1 - Ax) 2 + z 2 (1 - Ax) (1 - 12x) + 2(1 + Ax) } , (13) 



do gg 7is 3 VT^Ax~F 2 (n) (tt 2 + AT 2 (n)) 



dcosO 1024 

x {3z 4 (1 - Ax) 2 + 2z 2 (l- Ax) (5 + 12x) + (1 + I2x) (3 + Ax) } , (14) 



F(n) = s^/M^ 2 , T(n) = I(M S /V§), (15) 
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where x = m 2 z /s, z = cos9, a = e/(4 sin9w cos9w), and b = — 1 + 8sin 2 9w/3 or —1 + 
4 sin 2 9w/3, for up-type and down- type quark, respectively Note that the factors F{n) and 
T(n) satisfy 

k 2 D(s) = 87rJ !r (n)(7r + 2iT(n)). (16) 

In the numerical calculations we use the CTEQ6L1 PDF set and choose the factoriza- 
tion scale to be the Z boson mass. To estimate the uncertainty of the tree-level calculations, 
we vary the factorization scale from m>z to 1 TeV, and find that the uncertainty of the total 
cross section is about < 50%, corresponding to a rather small shift of Ms- Fig. [2] shows the 
total cross sections in the LED model as functions of Ms- We can see that the cross sections 
drop quickly with the increasing Ms, but can still reach hundreds of fb for Ms = 4TeV. At 
low Ms values, the gg fusion channel is dominant, while for large Ms the qq annihilation 
contributions become important and are comparable to the ones of gg fusion channel. 



IV. MONTE CARLO SIMULATIONS 



We only consider the 4 charged leptons final state since at such high energy scale the two 
jets from highly boosted Z boson decay are hard to distinguish and the jet final states will 
suffer from large backgrounds also. The partonic level events are generated with COMPHEP 
4.4 17] and PYTHIA6.4 [3] is used to treat initial and final state radiations. We smear 
the energy and direction of the charged leptons by a Gaussian distribution with the errors 
given by [l9|] 



AEi/Ei = 0.1/ WGeV © 0.007, 
Arn/m = A^j/0, = 0.001. 



(17) 



Here we do not distinguish 
acceptance cuts as follows 



between electrons and muons for simplicity. We impose the basic 



19] 



p l T > 15GeV, \r)i\ < 2.4, and AR U > 0.1. 



(18) 



For the reconstruction of the two Z bosons we require the lepton combination of minimizing 
the following quantity 



A = a/ (m hl2 - m z ) 2 + (m hh - m z ) 2 . 



(19) 
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In Figs. [3] and |4j we show the 4 lepton invariant mass distributions and the leading 
lepton pt distributions, respectively. It can be seen that the LED signals prefer the high 
invariant mass and high p? regions. And for Ms = 1.5TeV, the LED contributions can 
change the shape of the distributions significantly, while for Mg = 5TeV, the changes of the 
distributions are relatively small, leading to long tails in the distributions. Fig. gives the 
velocity or boost factor (3 distributions of the rest frame of the Z pair. We can see that the 
SM distribution will have two peaks in the forward and backward region since the Z boson 
pair can only be produced through qq annihilation in the SM. But for the LED contributions, 
gg fusion is dominant at low Ms values, thus the distribution will be flat and drop in both 
ends. In Fig. we display the polar angle distributions of the reconstructed Z boson in the 
rest frame of the Z pair. Compared with the SM case, the LED signal falls down in both 
ends since the leptons from the decay of highly boosted Z boson in these regions fail the 
rapidity cut. 

In order to further investigate the LHC reach of Ms, we introduce an invariant mass cut 
m/a > 1.5TeV and use a simple event counting method to calculate the signal significance. 
The Z boson pair production cross section in the SM is normalized to 17pb based on the 
next-to-leading order result of Ref. [9|. Table dgives some results of the cut efficiencies, signal 
significance and total cross sections after all the cuts. It can be seen that the invariant mass 
cut reduces the SM production rate significantly while only has small influence on the LED 
signals, especially for large Ms values. In Fig. [7] we show the LED signal cross sections as 
functions of Ms after all the cuts, and also the LHC reach of Ms- Generally speaking, for 
the Z boson pair production process, the LHC can probe the Ms values up to 4.3 ~ 5.6 TeV 
for n = 6 ~ 2 at 3a level, assuming an integrated luminosity of 100 

It should be noted that since we cut off the subprocess energy at some scale, we didn't 
include the LED contributions beyond the cutoff in our calculations, because they de- 
pend on the unknown ultraviolet completion of the LED theory. For low Ms values, like 
Ms = 1.5TeV, these contributions may be comparable or even larger than the ones we have 
considered, which may change the total cross sections and all the distributions. On the other 
hand, for high Ms values, like Ms = 5TeV, the LED contributions beyond the cutoff are 
small due to the suppression of the parton distribution functions and the limited energy of 
the LHC, and can be neglected. So these contributions will not change the above results of 
LHC reach for Ms, which is about 5 TeV. 
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cut efficiency 


M s = 2TeV 


M s = 3TeV 


M s = 4TeV 


SM 


n=2 


n=3 


n=6 


n=2 


n=3 


n=6 


n=2 


n=3 


n=6 


basic cuts 


0.73 


0.73 


0.76 


0.76 


0.78 


0.79 


0.72 


0.74 


0.75 


0.35 


m Al > 1.5TeV 


0.25 


0.30 


0.42 


0.51 


0.59 


0.66 


0.49 


0.56 


0.61 


0.00027 


total cross section (fb) 


27.2 


18.1 


8.2 


6.49 


3.87 


1.36 


1.22 


0.69 


0.22 


0.046 


S/^S + B 


52 


42 


28 


25 


20 


12 


11 


8.1 


4.3 





TABLE I: Cut efficiencies, signal significance and total cross sections after all the cuts with different 
parameters, assuming an integrated luminosity of 100 fb^ 1 . 

V. CONCLUSIONS 

In conclusion, we have studied the Z boson pair production mediated by the KK graviton 
in large extra dimensions at the LHC. Since the scattering amplitudes grow quickly with the 
energy, we used the partial wave unitarity to discuss the constraints on the process energy 
scale. We found that the LED contributions can enhance the Z boson pair production cross 
sections significantly when the fundamental scale Ms of the large extra dimensions is up to 
several TeV, and the kinematic distributions of the LED signals are greatly different from 
the SM ones through Monte Carlo simulations. Finally, we investigated the LHC reach of 
Ms and found that the LHC can probe the M s values up to 4.3 ~ 5.6 TeV for n = 6 ~ 2 at 
3<T level assuming an integrated luminosity of 100 fb^ 1 . 
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FIG. 1: Partial wave unitarity constraints on the ratios yfs/Ms and A/ Ms, the allowed region is 
to the left of and below the curves. 



11 




0.01 -| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

1000 1500 2000 2500 3000 3500 4000 4500 

M s (GeV) 



FIG. 2: The total cross sections of the Z boson pair production from LED contributions as functions 
of M s at the LHC. 
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FIG. 3: The invariant mass distributions of the 4 leptons. 
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FIG. 4: The transverse momentum distributions of the leading lepton. 
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FIG. 5: The normalized velocity distributions of the rest frame of the Z pair. 



15 




FIG. 6: The normalized polar angle distributions of the reconstructed Z boson in the rest frame of 
the Z pair. 
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FIG. 7: The total cross sections of the LED signals after all the cuts as functions of Ms- The 
horizontal lines indicate the cross sections needed for a 3a detection of the signal. 
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